Small distance expansion for radiative heat transfer between curved objects 
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We develop a small distance expansion for the radiative heat transfer between gently curved 
objects, in terms of the ratio of distance to radius of curvature. A gradient expansion allows us to 
go beyond the lowest order proximity transfer approximation. The range of validity of such expansion 
depends on temperature as well as material properties. Generally, the expansion converges faster 
for the derivative of the transfer than for the transfer itself, which we use by introducing a near-field 
adjusted plot. For the case of a sphere and a plate, the logarithmic correction to the leading term 
has a very small prefactor for all materials investigated. 
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More than 40 years ago Van Hove and Polder used Ry- 
tov's fluctuational electrodynamics [ij to predict that the 
radiative heat transfer (HT) between objects separated 
by a vacuum gap can exceed the blackbody limit This 
is due to evanescent electromagnetic fields decaying expo- 
nentially into the vacuum. HT due to evanescent waves 
has also attracted a lot of interest due to its connection 
with scanning tunnelling microscopy, and scanning ther- 
mal microscopy under ultra- high vacuum conditions 0, 0] ■ 
The enhancement of HT in the near-field regime (gener- 
ally denoting separations small compared to the thermal 
wavelength, which is roughly 8 microns at room tempera- 
ture) has only recently been verified cxperimentally[5|,0]. 
Theoretically, HT has been considered for a limited num- 
ber of shapes: Parallel plates a dipole or sphere 
in front a plate (ol-flU , two dipoles or spheres [1, [iJ, [3 , 
and for a cone and a plate [Ij] . The scattering formalism 
has been successfully exploited [lol. ITEI - flsj in this context. 
Although powerful numerical techniques [13, 13 exist for 
arbitrary geometries, analytical computations are limited 
to planar, cylindrical and spherical cases [3. Alterna- 
tively, the HT between closely spaced curved objects can 
be computed by use of the proximity transfer approxima- 
tion (PTA) [13, [m, HO] , which exploits an approach that 
has been extensively used in the context of fluctuation 
induced forces (2]| (referred to as proximity force approx- 
imation): The HT between two parallel plates (per unit 
area), Hpp(S), for separation S is averaged over one of the 
(projected) curved surfaces. PTA is generally assumed 
to hold asymptotically for small separations, however no 
rigorous derivation appears available in the literature. 

Here we develop a gradient expansion for HT between 
closely spaced objects, which enables to rigorously jus- 
tify PTA and to quantify corrections to it in the limit 
of small separation. This method, which has been pro- 
posed for Casimir forces [2^-l2^. exploits the mapping of 
coefficients of a perturbative expansion on one side and a 
gradient expansion on the other. We carefully explore the 
limitations and subtleties of this method in the case of 
HT and propose a near-field adjusted plot that overcomes 
the possibly slow convergence of the expansion. 

Consider two non-magnetic objects with dielectric per- 
mittivities ei(w) and e2(w) at temperatures Ti and T2, 



respectively. The HT between them is symmetric and 
can be written as an integral over frequency uj, as 



H = 



dcu [n„(Ti) - n^(T2)] h, 



(1) 



where n^(T) = {e'"^/^BT _ jg Bose- Einstein 

weight. The spectral density h depends only on the clas- 
sical scatterin g am plitudes of the objects and their rela- 
tive positions ^M,- For the sake of simplicity, let object 
f be a semi-infinite, planar body, filling the space z < 0, 
whereas object 2 has a curved surface with the smooth 
height profile z = S'(x); x = {x,y) are Cartesian coor- 
dinates in the plane normal to the z axis. PTA is then 
calculated as 



HpTA 



(2) 



and analogously for the spectral density, where E is the 
projected surface. 

Different contributions [qJ to Hpp{S) are separately de- 
picted in Fig. [T] for two dielectric materials, SiC [25| and 
Si02 (computed using optical data [1^). Propagating 
waves give rise to an almost separation independent con- 
tribution, whereas the transfer due to evanescent waves 
depends strongly on separation. Both parts have in turn 
two distinct components from electric (E) and magnetic 
(M) modes. Notably, if S is the smallest scale, the 
evanescent i?-modc diverges as S~'^, while all other con- 
tributions approach finite values. For the two materials 
considered, the figure indicates that the contribution 
is dominant for 5' < 0.03//to and 5' < 0.3/im, respectively. 
In general, the range of this regime, sometimes referred 
to as the nonretarded regime, depends strongly on ma- 
terial properties [13]. Because the behavior of Hpp{S) 
is nontrivial in general, we restrict ourselves to the case 
Hpp{S) cx S-^, such that Eqs. ©-((n]) below are vahd. 
At the end, we give an outlook to larger separations. The 
expansion can equally well be performed for the transfer 
in Eq. ^ as for its spectral dendity, and we choose the 
latter due to its larger information content. The spectral 
density for two parallel plates reads accordingly 



(3) 
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FIG. 1. The evanescent electric magnetic (J/p 

and total propagating (Hpp) contributions to heat transfer 
between parallel plates per unit area, Hpp = Hpp'^ + Hpp'''^'' + 
Hpp, as a function of separations for Ti = OK and T2 — 
300K. Data are shown for two materials, SiC and Si02, and 
are normalized by the result for black bodies (where a — 



where the coefficient A^^ uniquely describes the material 
dependence of the HT in the considered range. Follow- 
ing Refs. [l^-H^i we assume (justified a posteriori) that h 
admits a local expansion in VS*, 



h[S] =hpTA[S] 



(fx hpp{S) \/S-\/S- 



(4) 



where /3uj is a geometry-independent expansion parameter 
and the dots represent higher order terms in VS". We 
will see below that Eq. (U) is an expansion in the ratio of 
separation to radius of curvature of the surface. 

The coefficient (3ujiS) is obtained from the perturba- 
tive expansion Q of h[S] for |s(x)| <C d, with s de- 
fined with respect to the point of closest proximity d, 
S'(x) = d + s(x): If both a perturbative expansion as 
well as the gradient expansion in Eq. (U) exist, the two 
expansions must have identical coefficients in their com- 
mon limit, i.e. the expansion of h in Eq. @ has to match 
the perturbative expansion at small momenta k, which 
is given by 



h[d + s(x)] = ao{d) + ai{d)s{0) 
d^k 



a2{k-d) |S(k)|^ 



(5) 



where s(k) is the Fourier transform of s(x) with in-plane 
wave- vector k. Thus, if it is possible to expand a2{k\d) 
to second order in k, we can match the expansion 

a2{k;d) = a^^\d)+a^^\d)k'^ + (6) 
to Eq. @ to get 

Pu.hpp{d)^af\d). (7) 



In the asymptotically close regime, (3^ is d-independent, 
as both hpp{d) and 0^2^^ (rf) are proportional to d~^. The 
perturbative calculation for h[S] is done by use of a plane- 
wave basis with in-plane wave-vector k (eventually iden- 
tical to k in Eq. ([5])) and modes p = E, M. The reflection 
matrix of the planar object 1 is diagonal in this basis and 
given by the Fresnel coefficients r^^-* (k). The correspond- 
ing coefficient R^^J, (k, k') for the curved object 2 depends 
on both incoming (primed) and scattered waves, and can 
be expanded in powers of s[28j (with k = |k|), as 

i?J,5(k,k') ^ (2^)2j(2)(k- kOv^i-'^fc) 

+ i?«(k,k')S(k-k') 



d^k" 

(2^)2 



<)(k,k';k")S(k-k")S(k"-k') 



PF 



(8) 



The expressions for B^^, and -BppV are given m 
Substituting Eg. (j8]) into the HT-exprcssion for corru- 
gated surfaces 1 15l| and expanding to second order in s(k) 
(2) 

yields a2 (d). Expanding the latter for small k yields 
a2{k;d) according to Eq. ([6]), and with Eq. ((T)) the de- 
sired f5c^ is determined. 

For the experimentally most relevant case of a sphere 
of radius R in front of a plate, Eq. (|4]) is evaluated to 
yield the following expansion in d/ R (where iS'(x) = d + 

i?(i - yr^wW)), 



Jii) 



Ref. 



h{d) 



2'kRK 



1 -(2/3.-1)^ log 



0(d"). (9) 



Here do is an unknown constant of integration. PTA in 
Eq. ^ corresponds to the result for /3w = (compare 
to Eq. dH)). Hence in Eq. ^ 2Ti\^\og{d/ do) represents 
a trivial correction to the leading 2TTRX^/d and appears 
due to the projection of the sphere's surface onto the 
xj/-plane. 

As the term of 0{d'^) in Eq. ^ is unknown from the 
expansion and might not be small compared to the other 
terms in Eq. ^ (we compare it to a logarithmic term), 
possibly leading to slow convergence in practical cases, a 
better quantity to consider is its derivative. 



1 + (2/3. 



, d 
'^R 



(10) 



which converges faster, as the constant term in Eq. ([9]) 
drops out. 

In order to test these results, we numerically compute 
the exact heat transfer between a sphere and a plate, 
by utilizing the scattering properties of sphere and plate 
in cylindrical basis [l4|. Implementing a multipole ex- 
pansion, a maximal order of Zmax = 2500 is used for the 
smallest separations d/R^ 0.004 shown to obtain precise 
data (including the derivative). We can hence validate 
Eq. ([TU|) (or Eq. ^) by finding /3„ in two independent 
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FIG. 2. (color online). The gradient expansion parameter 
as a function of frequency for SiC, computed via Eq. ((Tjl, 
(solid line), and as found by fitting exact data for a sphere 
in front of a plate to Eq. (|1U|I (data points). The dashed 
line shows hpp{d = 10nm)nuj{T — 300K), illustrating the 
dominant range for near-field transfer of SiC. The upper in- 
set provides three examples of fitting the exact data (data 
points) to Eq. (|10[l . where solid (dashed) lines omit (in- 
clude) the assumed next order term in Eq. (|10|l [29l |. The 
point styles used are in accordance with the corresponding 
frequencies in the main figure. In the considered regime 
(Eq. ((Sjl) these curves depend on R only through d/R. For 
uj = {0.5934, 0.5967, 0.6050, 0.6075, 0.6100}rad//im shown in 
the main plot, the values computed in the gradient expansion 
(/3„ = {0.774, 0.206, -0.086, -3.026, -2.156}) agree weU with 
the fitted ones ({0.773, 0.206, -0.105, -3.044, -2.112}). The 
lower inset shows Pu, for Si02. 



vifays. The main panel of Fig. [5] shows /J^^ computed by 
the gradient expansion for SiC, together with five values 
obtained by fitting of exact data to Eq. (fTO|) . We observe 
good agreement [29| . The plotted frequency region is the 
relevant interval for SiC, with one dominant resonance 
around uj = 0.6 rad//xm, as indicated by hpp{d = lOnm) 
multiplied by the weight 71lj{T = 300K) (dashed curve). 
The upper inset shows the result of Eq. ([TU)) using the fit- 
ted Pu: together with the exact data for three frequencies 
also depicted in the main plot. We have thus demon- 
strated the validity of the gradient expansion both quali- 
tatively and quantitatively. The lower inset shows for 
Si02 computed by the gradient expansion. 

The total transfer H follows from Eqs. ([T]) and 
where we still consider Hpp{S) oc S^^. For materials 
with jej not too large (e.g. Si02), this regime requires 
d <C {AtijAts} (At = hc/{kBT) is the thermal wave- 
length. At ~ 7.6/im at 300K. It sets the range of wave- 
lengths contributing to the transfer in Eq. ((T))). Integrat- 
ing Eq.(|n]) over frequencies lu, we obtain 



H{d) = 



2t:RX 



, d d 



+ 0{d°), (11) 



Figure [3] depicts the HT from a sphere to a plate as 
a function of d/R for Ti = 300K and T2 = 0. In or- 
der to eliminate the constant term in Eq. ljlip . the main 
panel depicts a near-field adjusted plot, where we show 
H{d) — H{d = 0.004i?). Such a plot takes advantage 
of the faster convergence observed for the derivative in 
Eq. ([TU|) , and allows to accurately describe the exact data 
by use of Eq. (jlip . and in the same manner, experimen- 
tal data could be compared to Eq. (fTTj) . We note that 
the reference separation (here d ~ 0.004i?) should be 
chosen in the range where the d/i^-expansion is valid. 
For both SiC and Si02, the exact data points agree re- 
markably well with Eq. pT|). over the whole range of 
separations depicted (up to d/R = 0.1). We also show 
PTA of Eq. which - again - is to the lowest two 
orders m d/R given by Eq. ([TT|) with ji — Q. Its devia- 
tion from the exact data is hence quantified by Eq. (fTTj) 
to 47r/3Alog[(i/0.004i?]. This is shown in the inset, with 
convincing agreement to the exact data [3l|. Based on 
Fig.[Tl we expect these curves to describe the HT for sep- 
arations up to d < 0.03/im and d < 0.3/im for SiC and 
Si02 respectively. E.g., for R = 1/im [R — 10/xm) for 
SiC (Si02), this corresponds to d/R < 0.03. 

The subleading term in Eq. (jlip is a logarithmic cor- 
rection. Surprisingly, this term is hardly noticeable due 
to the coincidence that /3 = 0.5119 and /? = 0.5241 for 
SiC and Si02, respectively (despite very different curves 
for for the two materials in Figl2]), and hence the 
prefactor (2/3 — 1) is very small (0.02 and 0.05). In other 
words, the second order term in Eq. (jlip is very small 
because the term with (3, coming from the second term 
in Eq. (|4]), almost cancels the one coming from the first 
term in Eq. . Due to this coincidence, which is special 
for the sphere-plate geometry, PTA predicts a logarith- 
mic term which in reality is almost absent. 

Other geometries can readily be analyzed with Eq. (|4|) . 
For example, for two spheres of radii Ri and R2, we find. 



Hid)^ 
-(2/? 
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log 
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d d_ 

Tr^ °^d^ 

0{d% 



(12) 



where, in contrast to the sphere-plate case, the logarith- 
mic term is pronounced for /3 0.5, as has been observed 
in numerical data (20| . Furthermore, for a cylinder of ra- 
dius R and length L in front of a plate, the expansion 
convergences better than Eq. ([TT|) . because the second 
term has no logarithm. 



H{d) 
L 



RX 



V2d3/2 



1 



2/3- 



0{d% (13) 



where /3 and A are computed from f3^, A^j , Ti and T2 [3C 



Geometries for which semi-analytical methods of scat- 
tering theory are not possible (e.g. paraboloids or ellip- 
soids), can also be studied with the aid of Eq. (|4]), and 
will be left for future work. 

For larger separations, the d-depcndence of the expan- 
sion parameter in Eq. ([7]) has to be taken into account. 
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FIG. 3. Near-field adjusted plot {H{d) - H{d = 0.004i?)) for 
a sphere in front of a plate for SiC (circles, left abscissa) and 
Si02 (squares, right abscissa), Ti = 300K and T2 = OK. Exact 
data, PTA and Eq. (Illf) are shown as data points, solid and 
dashed lines, respectively. The inset depicts the deviation of 
PTA from the exact data, where lines show the analytical form 
47r/3Alog[d/0.004JJ] from Eq. (fTTj) . In the regime considered 
(Eq.Q), all curves shown depend on R only through d/R. 



and Eq. ^ has to be evaluated numerically. This anal- 
ysis, which poses a number of additional challenges [s^ 
(as well as insights), is left for future work. 



To conclude, the heat transfer between smoothly 
curved objects can be systematically expanded for small 
separations. This expansion is well suited for the gradi- 
ent of the transfer, and a near-field adjusted plot allows 
accurate prediction of experimentally measurable data. 
If the separation is the smallest scale, the expansion can 
be performed analytically, where coefficients are evalu- 
ated numerically. 
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